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Abstract 

In the present work we demonstrate how to reahze Id-optical closed lattice experimentally, 
including a tunable boundary phase- twist. The latter may induce "persistent currents", visible by 
studing the atoms' momentum distribution. We show how important phenomena in Id-physics 
can be studied by physical realization of systems of trapped atoms in ring-shaped optical lattices. 
A mixture of bosonic and/or fermionic atoms can be loaded into the lattice, realizing a generic 
quantum system of many interacting particles. 
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Studies of one dimensional systems constitute an intense research activity both in exper- 
imental and theoretical physics. They are particularly interesting mainly because quantum 
effects are strongest at low dimensionality and peculiar phenomena emerge. Prominent ex- 
amples are the spin-charge separation in Luttinger liquids one dimensional persistent 
currents in mesoscopic rings ^, and transmutation of quantum statistics!^. Most of the 
approximate schemes working in higher dimensions break down in 1-d. Only for a restricted 
class of model Hamiltonians, physicalproperties can be obtained analytically resorting to 
powe.M .e.hn„ue. a. Be.Ke a^aj co,*..a, field Fo. more generic Id 

systems, numerical analysis is the standard route to extract physical information. Degen- 
erate atoms in optical lattices could constitute a further tool for the investigationsml, thus 

Q 

rediscovering Feynman's ideas[7| suggesting that an ideal system with a "quantum logic" 
can be used to study open problems in quantum physics. Precise knowledge of the model 
Hamiltonian, manipulation of its coupling constants, possibility of working with controllable 
disorder are some of the great advantages of atomic systems in optical lattices compared 
with solid state devices to experimentally realize Feynman's ideas. The upsurge of interest of 
the scientific community has been remarkable, and some perspectives disclosed by trapped- 
atom "labs" have been already explored: the observation of the superfluid-Mott insulator 
quantum phase transition^, the analysis of the To nks- Girardeau regime in strongly inter- 
acting bosonslSj, and the physical realization of a Id-chain of Josephson junctions 10| were 
relevant achievements for condensed matter physics. The two most widely used methods to 
trap and manipulate atoms are based on the conservative interaction of atoms with either 
magnetic fields or with far off-resonant laser beams. For our purposes the magnetic trapping 
potential has a parabolic symmetry. Laser light interacts with the atomic induced dipoles 
creating attractive or repulsive potentials depending on the sign of the detuning A from 
resonance This can be used to create different potentials for different atoms, but with 
a single tunable laser beam. Notice that no light absorption occurs in creating the potential; 
therefore the medium can be considered transparent to the laser. 

So far open optical lattices have been studied. This constitutes a limitation of opti- 
cal apparata since a variety of studies for finite Id lattices with Periodic Boundary Con- 
ditions (PBC) exists in the literature, that cannot be accesed with them. In the same 
way as Gaussian laser beams are useful to produce open optical lattices, we shall take 
advantage of the rotational symmetry of Laguerre-Gauss (LG) laser modes to produce 
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FIG. 1: The optical potential resulting from the interference of a plane wave with an LG mode 
with L = 14, p = 0. For p ^ the potential is virtually unaltered 

closed optical lattices. LG beams, obtained experimentally making use of computer gen- 



erated holograms 



12| |. have already been used in the field of ultra-cold atoms[12i]. A LG 



mode with frequency u, wave-vector k and amplitude Eq propagating along the z axis 
can be written in cylindrical coordinates {r,ip,z) as E{r,ip) = EQfpi{r)e^^'''e^^'^^~^^\ 

fpii^) = ^^'|/hi ^'''"^p (^^) ^ ^ V^r/rQ, where is the waist of the beam, 

and Lp' are associate Laguerre polynomials L^{x) := {—)"^d'^/dx"^[Ln+m{x)], Ln+m{x) being 



the Laguerre polynomials themselves. The numbers p and / label the radial and azimuthal 
quantum-coordinates, respectively. The lattice modulation is obtained by interference of a 
LG beam with a plane wave Eqc'^^^^'''^^: in the far field, the interferogram is periodic in (p 
with I wells. For even / a perfect Irf-ring with L = I lattice sites is obtained. By reflecting 
the combined beam (LG beam plus plane wave) back on itself one achieves confinement also 
along z. Indeed a series of disk shaped traps are obtained. We point out that tunneling 
between the disks can be made much weaker than the corresponding tunneling within 
each ring adjusting vq/X (z. e. focusing the LG beam). Such a parameter depends 
monotonously only on L; for L > 15, t^/t^ < 10^^ can be achieved with Tq/A ~ 100. The 
resulting lattice potential (see Fig.((T))) is described by 

Viatt = 4^0 [1 + /pi + cos(/y.)] cos {kzf (1) 

Note that, contrary to what was done in here we need the laser frequency to be tuned 
below the atomic resonance since we want to trap atoms into the ring. For example with a 
laser intensity of / = 5W/cm'^ and A = —lO^MHz the potential wells would be separated 
by a barrier of ~ 5fiK much larger than the chemical potential of a standard condensate 
(whose temperature can reach few nK); with these parameters the scattering rate is -C 



1 photon/sec. It is worth noting that, because of the relation: Lm~™'(r) ^ Hn[{x — 
y) / \^]Hm[{x + y) / ^/2], LG modes can be reahzed also from Hermite-Gauss modes (modulo 
a 7r/2 phase change). Such a "mode-converter", realized experimentally in ij], can switch 
from an open to a closed lattice potential with the same periodicity and L. As we shall 
discuss further, this device might be useful in the experiments. 

We have just illustrated how to realize an optical lattice with PBC. Now we show how to 
twist them. The task can be achieved by applying an external, cone-shaped magnetic field 
B = B^pB^ + BzBz. In this way the atomic magnetic dipoles fimp experience a field varying 
along the ring, eventually equipping the periodic lattice by a twist factor: \I' e*'^™f'\I' at 
each „M„g, * bemg a genenc wave function. The pha=e facto. = ™,.cos., „,th 
tan^^ = B^p/Bz, is the analog of the Berry phasejly] of the two state system corresponding 
to the Zeeman splitting of the hyperfine atomic ground states; the role of time is played by 
the angle (p. We can adjust (pmp using an additional laser beam (with a suitable frequency), 
relying on the AG-Stark shift: AEifriF), where the function Ae depends on the intensity of 
the laser and on the Glebsch-Gordan coefficients corresponding to the matrix element of the 
electric dipole interaction energyj^]. The resulting phase twist is $0- = AEimp) +mF'n- cos 9 
where a = mp. Whereas boundary twists induced by a magnetic field pierceing the ring are 
"symmetric", = our protocol realizes IcZ-models with a tunable $0-; thus opening 
the way to novel investigations discussed below. 

For ORG, $0- can be "gauged away" completely from the system. In contrast, the bound- 
ary phase cannot be eliminated for closed loops and alters the phase diagram of the sys- 
tem [itI. Infact $o- emerges from the sum of site dependent phases causing an increase 
of the velocity field (oc to the tight binding amplitude t) that, in absence of dissipation, 
may set a persistent current. Therefore different regions in the phase diagram are identi- 
fied depending on the dynamical response of the system by perturbing The effect is 
reflected in the curvature of the A/'-particle energy levels En respect to the phase twist: 
Per = L"^ YlnPn \En{^a) " -E'n(O)] /{Ut^D, whcrc j9„ = c'^^"" /Z are the Boltzmann weights. 
For (spinless) bosons p+ = p_ = p is proportional to the superfluid fraction. Persistent 
currents are studied analyzing the charge stiffness P+ + P~ (for electrons, it is the 

zero frequency conductivity or Drude weight); a non vanishing Df. sets a persistent current, 
visible by releasing the condensate for a time much longer than the typical atomic oscillation 
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FIG. 2: Interference pattern for condensates released by the lattice obtained resorting to the 
analog of light diffraction from a circular grating[l^. The figures show the square of the order 
parameter \'4>{kx-,ky)\'^ = |?/;o(A;i:, fcj^) ^^Cq^ cos [«(A;j; cos(27rjyL) + /c^ sin(27rj/L) + </)j)] p. On the 
left = for all the condensates; on the right the atoms move along the ring with velocity oc V(/>; 
the interference pattern reflects a loss of matter at the trap center caused by centrifugal effects. 

period in the lattice wells. Then the spatial distribution of the condensates |\E'(r = kt,t)| 
is indicative of the initial atomic momentum distribution |\t'(k, 0)|[18|; in particular the 
phase difference between atoms trapped in different sites, produces characteristic interfer- 
ence patterns in the released condensates. In Fig. Q we show such a pattern for condensates 
released from the potential of Fig. in mean- field approximation (see also Fig.Q). Su- 
percurrent/superfluid fractions can be studied looking at the response of the system under 
imprinting of a dynamical phase (j)5t to the atomic wave functions, flashing the atoms 
with an additional Gaussian laser beam (can be much closer to resonance than those creating 
the potential) with a waist larger than the LG mode and with (y9-dependent intensity. The 
time 5t must be too short to induce atomic motion by absorption during the pulse). 

The case = — $_ is useful to study the spin stiffness Dg oc p+ — p_ indicating long 
range spin correlations in the system (for charged particles Ds would be proportional to the 
inverse bulk spin susceptibility^]). Generic values of $+ 7^ can be seen also as a 

result of certain correlated-hopping processes (on the untwisted models) [3] and correspond 
to more exotic cases that, as far as we know, have not been realized yet in physical systems. 
To be specific we consider A/" fermions described by the H ubbard model with particle-density 
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FIG. 3: The zero temp, momentum distribution for fermions with Hubbard dynamics is presented: 
\'^{k,,ky)\^ oc \w{k^,ky)\^j:i,je'^'^'''''''^Ek^ e*'^^(*»-*^)(nfc^); M/L = 32/16 {M/L = 1 is the less 
advantageous case to discern the effects of <I> at finite size, since the metal-insulator transition 

n 

strongly suppresses -Dc[2a|)j ('T-fc^) is calculated perturbatively, at second order in U/t. For <I> = 
(left), k^ = {-7r(AA- . . . 7r(AA - 1)/L}. For $ / 0, A;^ = {-7r(AA - + . . . 7r(AA - 
l)/L + ^/L}; the asymmetry in \"^{kx, ky)\'^ is due to the offset of {nk^) caused by ^. 

modulated kinetic energy 

i, (a) = t exp [i7,(a) + i J](a,-z(a)Arz,_, + A,- z(a)iVz,,) j , (3) 

I 

where Cj^^^s are fermionic operators, and Ni^^ := cj^Ci^^. U = Ti'bg f c/x|w(x)|^/m and 
t = J dxw{x)[—^'V^ + Viattlw^x + a), {bs, a, and w{x) indicates the scattering length, 
the lattice spacing, and Wannier functions respectively) play the role of the Coulomb and 



hopping amplitudes respectively; /ij^o- is of the order of the Bloch band separation j2l|; the site 
dependence can be achieved by tuning the magnetic confinement out of the symmetry axis 
of the optical ring. For the model (jSj) in a closed lattice, (jH]) can be gauged away everywhere 
but at the boundary; therefore (0), Q is equivalent to the ordinary Hubbard model, but 
with twisted BC0. The phase twist is <l>^ := 0(a) + 05-^I(a)7V_^ + 0V^|((t)(7V; - 1), where 



^) = E,=iai,m(^) . = E,=i (7j((^) + 0++(^) = E ^^=1 Aj^M + 

jT^m — 1,771 



^m,m-i(cr) + ^m-i,m+i(cr). Heuce, loading the Hubbard model into the twisted ring effec- 
tively leads to the physical realization of the model ©• To point out the effects of U 
(smearing of the Fermi distribution with algebraic singularity at kp) in the persistent cur- 
rent, |\l/(k)p is calculated for the Hubbard ground state at small U/t, and with $o- = 0(o"), 
= = (see Fig. ©). 
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The proposed setups could be used to study several issues in one dimensional systems. 

I. The concept of conformal invariance plays a central role in 1 + 1 dimensional critical 
phenomena: universality is characterized by a single parameter, the conformal anomaly c. 
The physical meaning of c resides in the concept of Casimir energy, namely the variation of 
the vacuum energy density to a change in the BC. For PBC it was shown [5J that the finite 
size correction to the bulk ground state energy is related to c: SpBc — £buik = —ncv/GL; 
resorting the modular invariance this correction should be visible in the specific heat of the 
system, at low temperature: C(T) = ttcLT/3v, for each collective mode of the system; the 
speed of sound v can be extracted from the dispersion curve, at small k: v = AS/Ak, for 



sufficiently large L (for the XXZ model, numerical analysis suggests that L > 15[25|). Except 
for integrable models, it is hard to measure or even have numerical estimates of c in solid 
.a.e sy.te,n|.H. W.tK *e presented setup, fo. MgK.v co„..ol.ab.e loaded .ode. these 
measurements can be done with unique accuracy. Both C(T) and AS/Ak can be measured 
following the techniques employed by Cornell et al. . To discern finite size effects in C (T) 
the PBC to OBC converter, discussed above, could be a valid tool. Indeed, the finite size 

co„.ect,o„ to ... fo.. OBC is also p.opo.iona. to c, but w.tK a Offerer, eoeffleieutfl . Then: 

8L 

cv = — {£pBC — £oBc) + where JF5 is the bulk limit of the surface energy that, being 

TT 

non-universal, can be fixed by performing the measurements for different L. (mimicking 
a "finite size scaling analysis"). Remarkably, both the energies £pBC and Eqec might be 
accessible measuring the second moment of the velocity of the released condensate (28l|. 
II. A general model we can engineer in the ring shaped lattice is 

H = Hbh + Huub + Hj (4) 



where Hbh is the Bose- Hubbard Hamiltonian 
boson interaction 



and Hi describes a density-density, fermion- 



By tuning A within the fine structure of the fermionic atoms, a spin 
dependence can be inserted in the hopping amplitude of the Hubbard model: t t„. At 
M /L = 1 and t„ <^ U the Hubbard ring effectively accounts for the physical realization 
of the twisted XXZ model with anisotropy 7 = (t^ + t^)/(2t+t_) and external field h = 
4 crt^//io-l21]. Loading quantum systems described by Hamiltonians of the type Q in 
lattices with twisted BC could serve to study charge and spin stiffness in physical systems 
with tunable interaction and/or disorder. For example a mixture of ^^Rb and "^^K atoms 
constitutes an ideal system to check the recent experimental evidence suggesting that the 
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supersolid order [s^ would be effectively favoured by the insertion of fermionic degrees of 
freedom into homogenous bosonic systems. The off-diagonal long range order manifests in 
superfiuid currents. Jumps between non-vanishing supercurrents should reveal the existence 
of the supersolid phase 3]| . This should be accompanied by a macroscopic occupation in the 
condensate at a non vanishing wave vector (~ 7r/(na), n > 2) signalling the charge-density- 
wave instability js^. The two condensates should be traced in the interference fringes. 
It was proved that exactly solvable twisted Hubbard/XXZ rings 0, 33 1 are equivalent 



to untwisted models for particles with intermediate statistics; this results in modifications of 
the exponents of the (low energy) correlation functions [ifll. The spatial profile of the latter 
might be detected by photoassociation techniques, as suggested in [3^ . 

Ill Another interesting issue we can study is the conjecture 3^ that Poisson or Wigner- 
Dyson level-statistics manifest in that the thermal Drude weights have qualitatively dif- 
ferent slopes for integrable (smooth algebraic temperature-decrease, universal behaviour of 
D(T)/ D{0)) or non-integrable (sharp, non-universal suppression of D(T)/D{0)) systems. 
Due to the precise knowledge of the model-Hamiltonian under analysis we can address the 
problem directly in a physical system. For example, we could consider "^^K pme-X X Z rings 
with twisted BC, for different L's; using the Feshbach resonance one could tune bg ~ 2a; 
the resulting XXZ model with next-nearest neighbor density- density interaction is non- 
integrable (another way is to destroy the integrability introducing disorder into the ring by 
site-dependent hj). In short: integrability can be switched on and off by tuning the Feshbach 
resonance (or adjusting the energy offsets hj). The presence of persistent currents can be 
detected along the lines described above (see Fig.Q, (jS)))- Numerical investigations for the 
XXZ model suggest that the effect should be visible for T/L > O.I7 j^]. 

In summary we have suggested a number of protocols to realize closed rings of many 
quantum particles, by optical means. This is possible by employing slight variations and 
combinations of techniques already developed within the current experimental activity in 
atomic physics. We have discussed how several open problems in condensed matter physics 
can be enlightened by such a setup. We finally observe that the clockwise/anticlockwise 
currents in a few-wells-ring constitutes a controllable two state-system analogous to the 
fiux-qubit realized by a SQUID. As the current is neutral, the corresponding decoherence- 
rate is much lower compared to solid state devices (charged currents). Information transfer 
could be mediated by an induced-dipole-dipole atomic interaction. 
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